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ABSTRACT 

A search for O VI 1031.93, 1037.62 A at redshifts corresponding to Lya lines in the z em ~ 2.4 QSOs HE1 122- 
1648 and HE2217-2818 reveals that a substantial fraction of those with H I column densities logA^(HI) > 14.0 
(cm -2 ) are highly ionized and show some heavy element enrichment. If these two sight lines are typical, then 
the O VI systems contain a cosmologically significant fraction of the baryons and the metals in the universe. For 
most systems the temperatures derived from the line widths are too low for collisional ionization to be responsible 
for the O VI lines. Photoionization models with a substantial hard ultraviolet flux can reproduce the observations 
for densities that are in good agreement with a model assuming local, hydrostatic equilibrium and heavy element 
abundances in the range ~ 10~ 3 -10~ 2 solar. Photoionization by a UV flux much softer than that predicted by Haardt 
& Madau (1996) for a background dominated by quasars can be ruled out. Finally, we find one system with a very 
low H I column density for which both photoionization and collisional ionization models yield a metallicity close 
to solar and a density that is inconsistent with gravitational confinement, unless the gas fraction is negligible. 

Subject headings: cosmology: miscellaneous — galaxies: formation — intergalactic medium — quasars: 
absorption lines — quasars: individual (HE1 122- 1648, HE2217-2818) 



1. INTRODUCTION 

There have been a number of searches for heavy elements in 
the low density intergalactic medium over a range of redshifts, 
motivated in part by the need to investigate possible enrichment 
mechanisms and to investigate the spectrum of the integrated 
ionizing background. The intergalactic medium is observable 
through the presence of Lya and heavy metal absorption lines 
in QSO spectra at redshifts up to that of the background QSO. 
Cowie et al. (1995) showed that for many of the high H I col- 
umn density systems (Nm <; 10 14 5 cm" 2 ) C IV is also present. 
Cowie & Songaila (1998) and Ellison et al. (2000) extended 
this analysis by a pixel optical depth technique to find that a 
similar C IV/H I ratio applies at lower column densities at red- 
shifts z ~ 3. The redshift dependence o f the Si IV/C IV column 
density ratio Songaila & Cowie 1996 suggests that there is a 
change in the nature of the ionizing background at z ~ 3, though 
a further in vestigation provide s little evidence for any redshift 
dependence Boksenberg 1997. 

While these early investigations concentrated on 
C IV 1548.20, 1550.78 since it has the strongest lines falling in 
the spectral region clear of contamination from the stronger 
Lyman forest absorption, at higher redshift photoionized 
OVI 1031.93, 1037.62 is expected to be a much better probe of 
the heavy element enrichment in the low density intergalactic 
medium (Chaffee et al., 1986; Rauch, Haehnelt & Steinmetz, 
1997). However, at high redshifts the presence of Lyman for- 
est lines makes the detection of the O VI lines very difficult. 
A thorough search by Dave et al. (1998) at redshifts ;> 3 did 
not reveal any convincing candidates. However, Schaye et al. 
(2000) have shown using a pixel optical depth technique that 
O VI is generally present in the Lyman forest absorption sys- 
tems seen in quasar spectra at redshifts 3. O VI absorption 
systems at low redshifts (z <J 1) are common (Buries & Tytler 



1996; Tripp, Savage & Jenkins, 2000), and a small number of 
such identifications have been found towards the z= 1 .73 QSO 
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At redshifts z J> 2 the O VI lines are in a wavelength region 
where they are measurable using ground-based facilities, and 
at z ~ 2 the Lyman forest confusion is less than at higher red- 
shifts. Therefore, it should be easier to detect the O VI lines at 
z ~ 2 than at higher redshifts. Here we describe the results of 
a search for identifiable lines from the O VI 1031.93, 1037.62 
doublet in specific absorption systems towards the z = 2.412 
QSO HE1122-1648 and the z = 2.406 QSO HE2217-2818. 
We investigate heavy element abundances and the nature of 
the ionizing mechanism for some representative cases where 
the H I, C IV and O VI lines could arise in the same regions. 
There are also systems similar to those reported by Simcoe et 
al. (2002), where there are clear velocity differences between 
these ions. 

2. OBSERVATIONS AND ABSORPTION LINE FITS 

The data reduction and analysis for HE22 17-2818 have been 
described by Kim, Cristiani & D'Odorico (2001), who also 
show part of the spectrum. Some spectral regions and results 
for HE1 122- 1648 are described by Kim et al. (2002). The 
spectral range covered by the analysis is from 3050 A to the 
C IV emission line in both cases. 

The continuum was initially determined as described by Kim 
et al. (2002), and the absorption lines fitted by multiple Voigt 
profiles, using the VPFIT program described by Webb (1987), 
Rau ch et al. (1992) as updated by Carswell et al . 
(see http://www.ast.cam.ac.uk/~rfc/vpfit.html). The whole of 
the available spectrum was fitted so as to provide a self- 
consistent set of line profile parameters for each ion. Except 
where noted otherwise below, different ions with components 

1 The data used in this study are based on public data from UVES Commissioning and Science verification at the VLT/Kueyen telescope, ESO, Paranal, Chile, 
institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK 

3 School of Natural Sciences, Institute for Advanced Study, Einstein Drive, Princeton NJ 08540 
4 European Southern Observatory, Karl-Schwarzschild-StraBe 2, D-85748 Garching bei Munchen, Germany 
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in the same system were not constrained to have the same red- 
shift, and their Doppler parameters were fitted independently. 
The main exceptions are the systems containing O VI described 
below where ionization levels were investigated. Each ion was 
fitted using all observed transitions which provided useful in- 
formation. 

Heavy element lines at wavelengths above the Lya emission 
line were fitted first, including components in the Lya forest 
region where necessary, and were used as indicators of where 
other ions might appear in the forest. Since we are interested 
in O VI 1031.93, 1037.62 over a range of redshifts, we have to 
be aware not only of the possibility of contamination by Lya 
absorption at lower redshifts, but also by Ly/3, 7 etc. Further- 
more, these contaminating lines affect not only the O VI lines 
themselves, but also the continuum around them. The strat- 
egy we adopted to allow for this, was to fit the Lyman lines in 
decreasing redshift order, taking care to include all previously 
fitted higher order Lyman lines when moving down in redshift. 
This procedure is generally satisfactory, but in a few cases a pu- 
tative low redshift Lya introduced to obtain a satisfactory fit to 
a higher order Lyman line at some redshift turns out to be pre- 
dominantly a different higher order Lyman line from another 
redshift system. In such cases the data were refitted with the 
misidentified Lya removed or modified. The net result after 
this iterative process is a set of Voigt profiles which fits the en- 
tire Lyman forest spectrum. Attempting to measure the O VI 
lines without doing this, and effectively using only O VI and 
Lya, would give misleading results in many cases. 

Fitting different transitions, particularly ones which are 
widely spaced in wavelength, revealed that in many cases the 
adopted continuum levels were not precisely correct, and local 
adjustments were needed to obtain a self-consistent fit. These 
continuum adjustments were made as part of the minimum \ 2 
parameter optimization within the VPFIT program. With few 
exceptions these continuum changes were less than 1 % of the 
originally adopted value, and tended to be correlated between 
neighboring fitting regions used. In light of this, a new con- 
tinuum was produced by interpolating over the continuum de- 
viations to provide a smooth, slowly-varying curve, which was 
then used to adjust the original continuum. The whole spec- 
trum was then refitted without allowing the continuum to vary 
to produce the final set of Voigt profile parameters. 

A sample list of the line fit parameters obtained in this way 
for HE1 122- 1648 is given in Table [j]. The table gives the ion, 
redshift z, Doppler parameter b in km s" 1 and log column den- 
sity N (cm -2 ) for each system found. The error estimates given 
against each quantity are from the diagonal terms of the co- 
variance matrix from the fitting procedure, and correspond to 
ltr estimates if the variables are uncorrelated. Where no error 
estimate is given, the variable involved has been tied to others 
on the same system. For the redshift z = 2.101764 example in 
the table, the redshifts of H I, C IV and O VI have been con- 
strained to be the same, and the Doppler parameters have been 
constrained so that b 2 = b^ h +bfa &lm , where b tur \, is identical for 
all ions and b 2 ha . m = 2kT/m, where m is the ion mass and the 
temperature T is constrained to be identical for all ions. The 
quantity b tm b is not necessarily true turbulence, but bulk mo- 
tion, such as the Hubble expansion of the absorbing regions, 
may be approximated in this way. We choose not to give ob- 
served wavelengths for the lines, since more than one transition 
from an ion was usually included in determining these param- 
eters. The full line lists for this object and HE2217-2818 are 



available as part of the electronic journal article. 

3. O VI SYSTEMS 

The fitted parameters for the systems where high ioniza- 
tion species may be present are given in Tables || and ||. For 
HE1 122-1648 the redshift range searched is 2.000 < z < 2.365 
and for HE2217-2818 2.000 < z < 2.360. At redshifts z < 2.0 
the signal-to-noise ratio in the regions of the spectra containing 
the O VI doublet is too low to allow meaningful measurements. 
We have imposed a high redshift cutoff at a redshift 4000 km s" 1 
below the emission redshift to avoid the region near the back- 
ground QSO. 

There are a number of ways in which absorption lines from 
different ions, which are thought to arise in the same system, 
can be fitted with Voigt profiles: the redshifts can either be con- 
strained to be the same or they can be allowed to vary indepen- 
dently, and the Doppler parameters can either be determined 
independently or be constrained to be consistent with a picture 
where the line is broadened by a mixture of turbulent and ther- 
mal motions. Since we are interested in ionization levels in 
single regions, we have chosen to constrain the redshifts to be 
identical and the Doppler parameters to consist of a turbulent 
and a thermal component as described in the previous section. 
Note that we implicitly assumed that all bulk/turbulent motions 
have a Gaussian velocity distribution. 

The la error estimates for the temperatures listed in Tables 
||and I are determined using an approximation based on small 
standard deviation normal distributions for the Doppler param- 
eters, so in some cases negative temperatures appear to be in 
the allowed range when turbulent velocities dominate. In the 
tables the redshift is left blank if it is constrained to be the same 
as that of the ion above it in the list. The Doppler parameters 
are given for all detected ions of each redshift component, with 
individual errors only for one ion as a guide, even though all 
line widths were included in determining these. 

The O VI 1031.93, 1037.62 lines are invariably in the Lya 
forest and are usually blended with higher order Lyman lines 
from higher redshift systems, so any putative identifications 
rely strongly on the presence of both lines with the correct dou- 
blet ratio. In many cases, as described individually below, the 
blending is dominated not by Lya, but by higher order Lyman 
lines or by some other species. In such cases the parameters of 
the contaminating lines are constrained by fits to other absorp- 
tion systems, and the measured strength of the O VI component 
is therefore more robust than if the contamination were dom- 
inated by Lya. Of course it is still likely that in some cases 
the O VI identification may be incorrect, with a chance super- 
position of Lyman lines and/or noise giving the appearance of 
the O VI doublet. Against this, there is a bias toward finding 
components with Doppler parameters less than ~ 15 km s , 
since broader ones are more easily confused with Lyman for- 
est systems. Our approach has been to try O VI in all possible 
cases, and retain it as a reasonably firm identification when a 
physically self-consistent fit can be obtained. The search for 
O VI was performed without regard to the strength of the cor- 
responding Lya line, since the neutral hydrogen fraction could 
be small. 

We can obtain an estimate of the number of spurious O VI 
identifications by using wavelengths offset from the real ones. 
A fake doublet with rest wavelengths 1035.0 and 1041.0 was 
chosen arbitrarily, but close to the real O VI doublet wave- 
lengths so that the signal-to-noise ratio in the data is similar. 
A search for fake O VI corresponding to H I with column den- 
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sities > 10 14 cm" 2 in HE1 122- 1648 revealed that in 8 of the 
24 cases the fake ion could be present. In the other 16 cases at 
least one of the fake O VI lines was demonstrably absent. This 
suggests that up to about 1/3 of the O VI identifications could 
be spurious. However, we note that many of our O VI identi- 
fications were then found to have detectable C IV at the same 
redshift, which increases the likelihood that they are real. 

It is interesting that for all but 2 of the 8 systems with 
log/V(HI) > 14.5 in HE1 122-1648 the data are consistent with 
the presence of O VI, and it may even be weakly present in one 
of the exceptions. In contrast, the presence of the fake doublet 
is excluded for 5 of the 8 systems. 

3.1. Comments on individual systems 

Below we present an atlas of all the absorption systems in 
which O VI could be present, and discuss each system indi- 
vidually. Even though most of the O VI lines occur in blends, 
the fits are much better constrained than the figures suggest 
because a large fraction of the absorption arises in higher or- 
der Lyman lines for which the parameters are determined from 
other hydrogen transitions. The uncertain cases are those for 
which there are blends with Lya components close to the O VI 
transitions. The Lya and other better-constrained components 
blended with the O VI lines are shown in Figures [j] to [i~2] along 
with the profiles for the lines in the individual systems. 

3.1.1. HE2217-2818 

z = 2.075 (Fig.0): 

Ly/3 shows the two-component nature of the H I clearly. 
The O VI lines at z = 2.075446 are not strongly 
blended with Lyman forest lines, but at z = 2.074740 
both the O VI 1031.93 and 1037.62 lines are seen as 
sharp components blended with broader Ly/3 from sys- 
tems at higher redshift and nearby strong Lya lines. 
C IV 1548.20, 1550.78 are present in both components, 
though in the higher redshift case there is a significant 
difference between the C IV and O VI redshifts. 

z = 2.180 (Fig.|): 

This system is probably more complex than the compo- 
nent structure given in the table would suggest. It shows 
ions from C II to O VI apparently coexisting in one com- 
ponent. We have not tried for a physically self-consistent 
model or estimated temperatures because of the complex 
structure. For all components in the complex the O VI 
lines are free of any significant blending from other sys- 
tems. 

3.1.2. HE1 122-1648 

z = 2.007 (Fig.|): 

The C IV 1548.20, 1550.78 lines are clearly broader than 
Si IV 1393.76, 1402.77, and are best fitted as two compo- 
nents, one at the Si IV redshift, z = 2.007 109, and another 
agreeing with O VI, z = 2.007205. Both lines in the O VI 
doublet are in a wavelength region where the signal-to- 
noise ratio is low. O VI 1031.93 is in the wing of a weak 
Ly6* and is otherwise free of significant contamination. 
O VI 1037.62 is seen as a distinct sharp component be- 
tween a Lya line at z = 1 .566566 and a blend of higher 
order lines along with Lya at z = 1 .567140. 
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FIG. 1 .— The systems at z = 2.074740 and 2.075446 in HE2217-2818 show- 
ing the profile fits to the H I, C IV and O VI regions (dot-dashed lines) against 
the data on a velocity scale centered on the redshift given above the figure. For 
H I and O VI each of the spectral regions is normalized to a continuum value of 
0.9, and an integer offset has been applied to separate them. In this case C IV 
scale has been multiplied by a factor three, so the zero levels for the C IV lines 
are 2.7 units below their continuum values. The individual components, here 
at z = 2.074740, 2.075398 and 2.075446, are shown as dotted lines in the lower 
six spectral regions. The upper two spectra show the O VI line regions with 
the features constrained by the presence of other transitions shown as dashed 
lines, and the blended Ly« as dot-dot-dot-dash. This example provides one 
reasonably clear and one blended suggested O VI identification. 
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FIG. 2.— The systems at z. = 2.180 in HE2217-2818 showing the fitted H I, 
C IV, O VI and Si IV profiles. The^C IV region has not been rescaled, and 
otherwise the details are as for Fig. hi 
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z = 2.030 (Fig.|): 

This system is unusual in that it shows C IV and 
NV 1238.82, 1242.80, but only weak Lya. Both lines of 
the O VI doublet are noisy and blended, and a marginally 
satisfactory fit can be obtained without their presence. 
O VI 1031.93 is in the wing of a blend of Ly/3 at z = 
2.048614 and Ly 7 at z = 2.215271, and O VI 1037.62 is 
in a blend of several higher order Lyman lines. There are 
no Lya lines near either component. 

z = 2.033 (Fig.|): 

There is a noise spike which obliterates the O VI 
1037.62 line, so the O VI identification at z = 2.033133, 
and parameter estimates, are based solely on a clean 
O VI 1031.93 line and its redshift coincidence with H I. 
CIV 1548.20 is not seen. 

z = 2.064 (Fig.|): 

The Lya profile is complex, and three components are 
present in C IV 1548.20, 1550.78. O VI is present in the 
z = 2.065681 component. The 1031.93 line is blended 
with Ly/3 at z = 2.084275. The O VI 1037.62 line is in a 
blend containing Lya at z = 1.616531 and several higher 
order Lyman lines from higher redshift systems. 

z = 2.080 (Fig.@): 

If the H I and O VI line widths are internally consistent, 
and the apparent double structure in each of the O VI lines 
is real, there appear to be at least two components, at 
z = 2.080162 and 2.080527. In a two-component fit each 
has a low estimated temperature but large turbulent line 
widths. The dominant blend with the O VI 1031.93 line 
is LyS at z = 2.346958, with parameters well determined 
from other lines in the Lyman series. The O VI 1037.62 
components fall between Ly/3 at z = 2.115206 and at 
z = 2.116542, and the wings of these lines affect the 
apparent continuum level against which the O VI line 
is seen. The blue wing of Lya at z = 1.629738 with 
b = 25.8 and log/V(HI) = 12.75 also has a small effect 
the O VI 1037.62 line profiles. 

It is also possible to fit this system as a single compo- 
nent, with the O VI lines having a Doppler parameter of 
39 km s" 1 and H I with a lower b (= 29 km s" 1 ), though it 
is not clear what this means. In this case log/V(HI) = 15.1 
and log/V(OVI) = 13.9. 

z = 2.102 (Fig.|): 

The structure around the O VI 1031.93 feature is domi- 
nated by Ly6 atz=2. 369161 shortward and C IV 1548.20 
at z = 1 .067739 longward of the observed wavelength 
of the line. Lya is present at z = 1.632438, but with 
logTV(HI) =13.07 and b = 49.7 this does not strongly af- 
fect the O VI line. The O VI 1037.62 line is significantly 
blended with Lya at z = 1. 647796 (log/V(HI) = 13.70, 
b = 28.1), and shows only as an extended blue wing to 
a feature at 3218.8A. 

z = 2.206 (Fig.|): 

C IV is weak, and the 1548.20 line may have a further 
velocity component. O VI 1031.93 is in the wing of 
a Lya line at z = 1.721641, and 1037.62 with Lya at 
z= 1.736637. 

z = 2.215 (Fig. [l^): 

Two components are detected, one with O VI line 



Doppler parameter b = 7.3 km s _I at z = 2.215271 and 
a more doubtful one at z = 2.215494 with broad lines. 
O VI 1031.93 at z = 2.215271 appears as a narrow com- 
ponent blended with Lya at z = 1.729154 (b = 35.9, 
log/V(HI) = 12.93), and for 1037.62 the higher redshift 
component is blended with Lya from a system with 
z= 1.744575,^ = 25.0 and log/V(HI) = 12.31. 

Z= 2.339 (Fig. [llj): 

In this case the presence of O VI rests solely on a weak 
1031.93 line at z = 2.339260. O VI 1037.62 at this red- 
shift falls close to a strong C IV 1550.78 at z = 1 .234266. 
The corresponding wavelength of Lya is significantly 
offset from the main component, and the H I column 
density at z = 2.339260 is very uncertain. The strongest 
neighboring H I components are included in Table ^| 

z = 2.352 (Fig.[l2|): 

Two components are measured at z = 2.352160 and 
2.352695. At each redshift the O VI 1037.62 line is suf- 
ficiently far from the nearest contaminant, a weak Lya 
in each case, that blending does not significantly affect 
them. At the lower redshift the O VI 103 1 .93 line appears 
as an extended wing of C IV 1548.20 at z = 1.234266. 
From the corresponding C IV 1550.78 it is clear that 
this extended wing cannot be due to C IV, even if the 
blended O VI 1037.62 at z = 2.339257 is absent there. 
At z = 2.352695 O VI 1031.93 is blended with Lya at 
z= 1.846002 {b= 18.1, log /V(HI)= 13.00). 
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FIG. 3.— The systems at z = 2.007109 and 2.007205 in HE1122-1648 
showing the fitted Hi, C IV, O VI and Si IV profiles. Ihe Si IV region has 
been scaled by a factor five. Other details are as for Fig. N. 
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HE1 122-1 648 z=2. 030096 
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FIG. 4.— The system at z = 2.030096 in HE1 122-1648 showing the fitted 
profiles in the regions of the H I, C IV, N V and O VI absorption lines. The 
C IV and. N V regions have been scaled by a factor three. Other details are as 
for Fig. hi 
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FIG. 5.— The system at z = 2.033133 in HE1122-1648 showing the fitted 
H I and O VI profiles. The C IV region has been scaled by a factor five, so 
the zero level is 4.5 units below the continuum. Note_the noise spike at the 
position of O VI 1037.62. Other details are as for Fig. hi 



HE1 122-1648 z = 2. 064789 
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FIG. 6.— The systems at z = 2.064412, 2.064789 and 2.065681 inHEmt 
1648 H I, C IV (scaled x5) and O VI profiles. Details are as for Fig. [J 
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FIG. 7.— The systems at z = 2.080162 and 2.080527 in HE1 122-1648 
showing tha fitted H I, C IV (scaled X 5) and O VI profiles. Other details are 
as for Fig. hi. 
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HE 1 122-1 64S z = 2. 10 \~r4 
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FIG. 8.— The system at z = 2.101764 in HE1122-1648 showing the fitted 
H I, C IV (scaled x5) and O VI profiles. Details are as for Fig. hi 
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Fig. 10.— The systems at z = 2.215271 and 2.215494 in HE1122-1648 

f owing the H I, C IV (scaled X 8) and O VI profiles. Details are as for Fig. 
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FIG. 9.— The system at z = 2.206460 in HE1 122-1648-showing the H I, 
C IV (scaled x5) and O VI profiles. Details are as for Fig hi 
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FIG. 1 1 .— The system at z = 2.339257 in HE 1 122-1 648-showing the H I, 
C IV (scaled X 8) and O VI profiles. Details are as for Fig. Ill 
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HE1122-1648 z=2. 352695 
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FIG. 12.— The systems at z = 2.352170 and 2.352695 in HE1 122-164 
showing theH I, C IV (scaled X 5) and O VI profiles. Details are as for Fig. |l 



4. OVERALL TRENDS 

One of the features of the O VI search is that in nearly all 
cases it is the high H I column density systems which show 
these lines. If logV(HI) > 14.5 then O VI is normally clearly 
present. This may be the explanation for the most obvious 
difference between the two sight lines; towards HE 1122-1648 
there are 10 distinct redshift systems with z > 2 showing O VI in 
some component, while towards HE22 17-28 18 there are only 
two. The numbers are too small for this difference to be statisti- 
cally significant, but it does correspond to the relative numbers 
of the high H I column density systems (14 and 6 respectively 
with z > 2 and logJV(HI) > 14.5). 
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FIG. 13. — The O VI Doppler parameter 
rameter for H I. Other details are as for 



The Doppler parameters for HI and OVI are correlated, as 
one would expect from the way the data were analyzed. The 
turbulence dominated systems have biO Vl)~ b(H I), and for 
those where thermal broadening dominates b(0 VI)~ 4fc(H I). 



It appears from 13 that both types of system are present, and 
that either turbulence dominates strongly or thermal broaden- 
ing does. For thermal broadening at temperature of 3 x 10 4 K 
b(H I)= 22 km s" 1 and b(0 VI)= 5.5 km s _1 . The group of 
systems with low O VI Doppler parameters have temperatures 
close to this value. 

Otherwise, it is hard to find any correlations at all. The 
O VI/H I column density ratio does appear to have an inverse 



correlation with the H I column density (see Fig. 14), so one 
might conclude that, as expected, the ionization level is linked 
with the inferred density. However, the O VI systems have to be 
quite prominent to be detected deep in the Lyman forest region 
of the spectrum and so most are near the detection limit. The 
apparent correlation may therefore arise simply because O VI 
lines corresponding to absorption systems in the lower part of 
the diagram are too weak to be detected. 




log N(HI) 

FIG. 14. — The O VI/H I column density ratio vs that for H I column den- 
sity, with lcr error estimates. Where no error bars are shown, the error range 
is within the symbol size. The open circles are for values where the O VI 
identification is doubtful. 



There are some correlations worth seeking since they could 
provide clues to the nature of the systems. We caution that the 
number of data points is small, so the absence of a significant 
correlation does not mean that the underlying reasons for an ex- 
pected trend are incorrect. One possibility is the O VI Doppler 
parameter against the O VI column density (Fig. |l5|). Here a 
trend could arise in several ways, such as undetected velocity 
structure in higher column density systems, cooling rates linked 
to metallicity, collisional ionization linking high ionization col- 
umn densities and the temperature, or just a selection effect 
since broad weak lines tend to be lost in the noise. No correla- 
tion is seen. Another is the inferred temperature vs H I column 
density, which might be expected if the temperature depended 
strongly on density. There is no such correlation discernible 
with the available data (Fig. Rq). 
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FIG. 15. — The O VI Doppler parameter, b, in km s 1 vs O VI column 
density. Other details are as for Fig. h4 
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FIG. 16. — fJh£ inferred temperature vs H I column density. Other details 
are as for Fig. Il4 

5. MODELING TFfE O VI SYSTEMS 

We have run CLOUDY 5 (version 94.00) |Ferl and 1999| mod- 
els for various forms of the ionizing flux using four example 
systems from HE1 122 - 1648 where the C IV and O VI appear 
to have the same velocity structure. These were chosen to cover 
a range of H I column densities and O VI/H I ratios and they are 
not necessarily the ones with the most reliable or accurately de- 
termined parameters. The details are given in Table |j. In this 
table the row containing the redshift identifies the system being 
considered, and gives the measured parameters from Table pi 
along with a density estimate based on the H I column density 
using a prescription from Schaye (2001), as detailed below. 

Subsequent lines in the table give the results of models for 
which the heavy element abundances and hydrogen density 
were the model parameters varied to give the best fit to the ob- 
served H I, C IV (where measured) and O VI column densities. 
The models were chosen to have all heavy elements present 
with solar relative abundances, but the heavy element/hydrogen 
ratio was allowed to vary to obtain a best fit to the data, using 
the OPTIMIZE command in CLOUDY. The log of the best fit 
C/H value relative to solar is given in the [C/H] column in the 
table. The other parameter which was allowed to vary to obtain 

5 http : / /www . pa . uky . edu/~gary /cloudy / 



a best fit is the hydrogen number density, «h- Since there are 
at most two independent ion ratios constraining the model pa- 
rameters, and two model parameters, a satisfactory fit is always 
possible. If C IV is undetected or unmeasurable then there are 
insufficient constraints to allow an estimation of both the metal - 
licity and the hydrogen density. 

The ionizing backgrounds chosen for these calculations are 
the Haardt-Madau model for the ultraviolet background radia- 



tion at z = 2 from quasars Haardt and Madau 199q , and simpli- 
fied power la w models. We hav e generally chosen a power law 
slope of -1 .8 Zheng et al. 1997 as the best estimate for a quasar 
background ionizing spectrum, and in some cases a power law 
slope of -1 .5 has been used to see how sensitive the results are 
to this choice. The effects of He II absorption are dealt with in 
a simple way by inserting a break so that the flux just above 
4 Rydbergs is 10~ B x the flux just below, with the same power 
law on either side of the break. For B we chose the values 
(no break), 1, and 2. Note that for a spectral index a, a break B 
corresponds to a softness parameter S = Thi/Thch ~ 4 a+1 10 B , 
where Thi and 1^11 are the ionization rates of H I and He II 
respectively. The z = 2 Haardt-Madau model has a softness pa- 
rameter S w 92, intermediate between the B = and B = 1 power 
law models. The model designation in the table gives the power 
law and break strength. For example, model "P1.8B1" has a 
power law slope of -1.8, and a break of 1 (i.e., a factor 10 in 
flux) at the He II Lyman limit. If the break is large (3, and oc- 
casionally 2 for the steeper power laws) there is very little high 
energy ionizing radiation and O VI is effectively suppressed. In 
those cases the O VI column density is close to zero and a fit to 
the data cannot be obtained within the operational range of the 
CLOUDY program. Such cases are omitted from the tables. In 
§[T3| we will show that photoionization by spectra with S 3> 1 2 
results in absorber sizes that are inconsistent with the data. 

The flux normalization was chosen to be logF„ = -20.31 
(erg cm" 2 s -1 Hz" 1 ) at 1 Rydberg in all power law cases. This 
corresponds to log/,, = -21.41 (erg cm -2 s -1 Hz -1 sr -1 ) and 
the photoionization rate for neutral hydrogen, log Thi « - 1 2 .02, 
matches the Haardt-Madau value if the power law slope is -1 .8. 
If a different F v value is preferred, then the results in the tables 
still apply, now for hydrogen densities rescaled so that the ion- 
ization parameter remains constant. 

We may transform the H I column densities to densities if 
we assume ph otoionization and (local) hydrostatic equilibrium 
Schaye 2001. This transformation is a function of the ionizing 
background, but we can use it here to estimate the background 
if we believe the models, or test the models if we believe the 
background, or use it to estimate the oxygen/carbon ratio com- 
pared with the solar value. For temperatures of ~ 5 x 10 4 K, and 
Thi ~ 10~ 12 02 , and choosing his default value for the fraction 
of the mass in gas ( f g = 0.16), Schaye's equation (8) linking the 
H I column density Nhi with the hydrogen density «h becomes 



log« H ~ -log/Vin-13.8. 



(1) 



We have used these density estimates to provide "observed" 
densities in tables ^ and ^ for comparison with those given by 
the best fit models to the ion ratios. Assuming ilbh 2 = 0.02, the 
mean hydrogen number density at z = 2 is ~ 5 x 10~ 6 cm -3 , with 
a corresponding Nm ~ 6 x 10 12 , so in all the cases considered 
here the regions have a significant over-density. 
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5.1 . Constraints based on HI, C IV and O VI 

For a given background ionizing flux the errors in the hy- 
drogen density and heavy element abundances given in Table 
U depend on the precision with which the column densities are 
determined. Given the assumption that the relative heavy ele- 
ment abundances are solar, changing the heavy ion/H I ratio ef- 
fectively determines the heavy element abundances, while the 
C IV/O VI ratio determines the ionization parameter and so the 
hydrogen density. For example, for the z = 2.101764 system 
in HE1 122-1648 the error estimate of 0.17 in logA^i becomes 
an uncertainty of about 0.17 in the heavy element abundances. 
The errors of 0.02 in logA^civ and logAfovi translate to an er- 
ror of about 0.03 in log«H- Thus, for a given background flux 
the hydrogen density can be quite precisely determined with 
the model assumptions we have made. However, the formal 
errors in the column densities and other quantities given in Ta- 
bles | and | are based on the assumption that the Voigt com- 
ponent model is correct in all respects, including in particular 
the number of H I components blended with the O VI lines. It 
is quite possible in the z = 2.101764 case that the O VI col- 
umn density is lower if there is additional H I blended with the 
O VI lines. If we choose, for example, to reduce the O VI col- 
umn density arbitrarily by a factor of two so logAfovi becomes 
13.83 without changing the C IV column density, then, for the 
Haardt-Madau background flux, the the best fit hydrogen den- 
sity becomes logn H = -4.00 and [C/H]= -1 .99. 

Collisional ionization models were also attempted, and the 
results for these are given the model designation "Coll" in Ta- 
ble Q The value in the temperature column is then the best fit 
temperature for the model. In no cases are the temperatures 
inferred from the line widths high enough to produce any sig- 
nificant collisionally ionized O VI. 

A further possibility is that different absorption lines arise in 
different parts of the cloud model, so our assumption of veloc- 
ity agreement may not be valid. For the CLOUDY models we 
have checked this point by examining the ion fractions through 
the cloud model, and find that the ionization fractions of the 
species we are modeling, H I, O VI and C IV, vary by at most 
10% through the cloud. This does, however, not preclude the 
possibility that different ionization levels come from different 
regions in a system with substructure. This point needs further 
investigation, but the fact that there is general agreement be- 
tween the ion redshifts when fitted independently suggests this 
is not a very important concern. 

There is also no guarantee that the C/O ratio in these sys- 
tems is close to the solar value of log(C/0) = -0.37, as as- 
sumed for the models. Since the enrichment is likely to have 
occurred from massive stars, the C/O ratio may well be smaller 
than the solar value. With examples of oxygen-rich supernova 
remnants as a guide (see e.g. Blair et al., 2000) suggesting 
that the C/O ratio can be about an order of magnitude smaller, 
we have run power law cases with oxygen enhanced so that 
log(C/0) = -1.37 while all other heavy elements have solar 
abundance ratios with respect to carbon. These are the mod- 
els labeled "OR..." in Table ^j, and the number in the "metals" 
column is then [C/H]. The net effect is an increase in the hy- 
drogen density to give the change in ionization then required to 
reproduce the observed O VI/C IV ratio. 

5.2. Temperature comparisons 

We now have temperature and density parameters inferred 
from the data for each system, and the same quantities, together 



with the element abundances, are predicted by the models. Be- 
fore exploring the agreement, or otherwise, between the two, 
it is worth examining the uncertainties in the two quantities in- 
ferred from the data in a little more detail. 

An important first question is: could the O VI arise from col- 
lisional ionization? As can be seen from a comparison of the 
temperatures inferred from the Voigt profile fits (Tables 0and 3|) 
and those needed for collisionally ionized models (TableQ), the 
temperatures inferred from the line widths are generally much 
too low to allow this as a possibility. However, the formal errors 
in the temperature estimates are based on the assumption that 
the decomposition into velocity components as separate clouds, 
each with a single temperature, is correct, and we know from 
models of the evolution of structure that this is an oversimpli- 
fied picture. And even if it were not, the decomposition into 
Voigt profile components is not unique. Since the components 
are optically thin the CLOUDY models based on the the relative 
column densities will still be broadly applicable, provided there 
are no large relative abundance variations between the compo- 
nents of what we treat as a single entity. However, the tem- 
peratures based directly on the line widths could be seriously 
affected. This is particularly true where there is a large turbu- 
lent component in the fit to the line widths, or when the lines 
are significantly blended. 

Temperatures T ~ 2.3-2.5 x 10 5 K are required if collisional 
ionization dominates, so these temperatures were imposed on 
each system. Acceptable Voigt profile fits to the O VI lines 
were found in all four of the test cases in Table ^j, either be- 
cause the Doppler parameters for the heavy elements are in any 
case large enough to accommodate such temperatures or, in the 
exceptional case of the z = 2.206460 system in HE1 122-1648, 
because other lines blended with the O VI adjusted to accom- 
modate the constraint and C IV is too weak to provide much 
additional information. The results are given in Table 0. Where 
the column density value is blank the error column gives the log 
of the la upper limit. 

The common feature to all of these high temperature models 
is that the Lyman lines associated with the O VI are too broad 
to be significant contributors to the observed H I line profiles. 
Additional, cooler H I components which contain the bulk of 
the H I column density are invariably required because the H I 
column density associated with the C IV and O VI must be low. 
Strictly speaking this new system would require remodeling to 
estimate abundances and collisional temperature, and then we 
should redo the fit and so on until the process converges, but 
we have not done this. The temperature is effectively set by 
the C IV/O VI ratio, and this can change somewhat when the 
H I constraint on redshift and Doppler parameter is changed. 
However, the observed range is satisfied by a narrow range of 
temperatures T ~ 2.2-2.5 x 10 5 K, so the predicted line widths 
will not be very different. The reduced H I column density then 
requires a model with higher heavy element abundances in the 
collisionally ionized component. If we require that the O VI is 
collisionally ionized, then each system showing O VI must have 
associated with it a cooler system with higher H I column den- 
sity and no detected heavy elements. Note also that in some of 
the cases given in Table [3] the C IV and O VI Doppler parame- 
ters are significantly less than 15 km s" 1 , and so do not allow 
the high temperatures required. 

At first sight even some of the photoionized models have 
temperatures which are too high. For example, for the z = 
2.101764 system the best-fit temperature is 1.3 ± 1.9 x 10 4 
km s" 1 , while the modeled temperatures are in excess of 4 x 
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10 4 K, which appears to rule out those models at roughly the 
2a level. However a look at the independently determined 
Doppler parameters for each ion in that system gives a hint to 
what may be happening. They are all much larger than the ther- 
mal Doppler parameters, and so bulk motion dominates the line 
widths. The best-fit turbulent component has a Doppler width 
of 29.1 ± 1.2 km s -1 , while even the H I thermal Doppler pa- 
rameter is ~ 15 km s" 1 . If we demand that the temperature for 
this O VI absorber is 50,000 K, for example, then a satisfactory 
fit to the data is obtained by introducing two velocity compo- 
nents at this temperature. Therefore, the temperature constraint 
is not one which is a significant concern for this system. 

The pair of systems at z = 2.080162 and 2.080527 cannot be 
dealt with in this way since the H I lines are themselves too nar- 
row to accommodate two components at a much higher temper- 
ature than about 20,000 K. It is possible to add a third, dominant 
H I component without associated heavy elements and obtain a 
satisfactory fit. However, the CLOUDY photoionization model 
equilibrium temperatures may not be the ones we should be 
attempting to match. In the hydrodynamic simulations adia- 
batic cooling as the universe expands is still important at these 
column densities. This cooling process is not included in the 
CLOUDY models, so temperatures they give are too high. 

5.3. Density comparisons 

If the absorbers are photoionized and in approximate, (lo- 
cal) hydrostatic equilibrium, then the densities predicted by the 
CLOUDY models should roughly agree with those predicted 
by the density /H I column density relation [equation (Q)]. From 
Table [| we can see that generally the models with harder far 
ultraviolet flux are preferred, i.e., those with no or a small 
change in intensity at 4 Rydbergs, and those with a flatter power 
law spectrum. The Haardt-Madau (1996) background spectrum 
produces models with densities to within about a factor two of 
the Schaye (2001) estimate. The pure power law models with 
no break and slope —1.8 also give density estimates in reason- 
able agreement, and for power law models with a slope -1.5 a 
break at the He II edge of a factor less than 10 is preferred. Any 
strong break in a power law spectrum gives densities that are 
significantly lower than expected. 

We can place a lower limit on the density from the bulk mo- 
tion component, b tur \,, of the line widths. The velocity width 
across an absorber of characteristic size I due to the Hubble 
flow is 

b H ~ H(z)£/2 = H(z)N H /2n H = H(z)N m /2n m . (2) 

Using equation 6 from Schaye (2001) for the neutral/total hy- 
drogen fraction of a highly ionized, optically thin plasma gives 

where T4 is the temperature in units of 10 4 K, is the H I ion- 
ization rate in units of 10~ 12 s" 1 , and we assumed (z,£l m ,flA) = 
(2.2,0.3,0.7). 

For systems with Nhi ~ 10 15 cm" 2 and densities less than 
~ 10 -4 5 cm -3 , the velocities involved are ~ 100 kms" 1 or more. 
Peculiar velocities can change this somewhat, but the Hubble 
velocity width is so sensitive to density that the upper limit 
£>turb ~ 40 kms" 1 (Tables |2] and effectively precludes the soft 
spectra models (..Bl and ..B2) given in Table ||. 

Note that both the Schaye (2001) and CLOUDY model den- 
sities have the same dependence on the normalization of the 



ionizing flux, so the agreement between the two does not con- 
strain this quantity in any way. The comparison is, as we have 
seen, sensitive to the shape of the ionizing background. So what 
we have shown is that, to obtain a consistent picture, the ioniz- 
ing background at z ~ 2 could have, from the choices we have 
explored, the Haardt-Madau form or a power law with a small 
jump at the He II edge. 

The density estimates of both the photoionization models 
and Schaye (2001) assume that the absorbers are illuminated 
only by the mean background radiation. If local sources of 
UV photons are significant, then the density will be underes- 
timated. The 2.007 system in HE1 122 -1648 illustrates the 
possible range well. The O VI is in a noisy region, so it was 
not selected for detailed analysis, but the presence of O VI in 
one system and Si IV in another with comparable H I column 
density suggests that the densities in the two components must 
be very different if the ionization source is the same for both. 
Table g shows that for power law ionization the density differ- 
ence is at least 1.5 orders of magnitude and large abundance 
differences are also required. Attempts to fit this complex as a 
single cloud fail. Unless there are large abundance anomalies it 
is not possible for O VI and Si IV to coexist in a cloud with H I 
column density logA^(HI) ~ 15.5. 

5.4. The high ionization system at z = 2.030096 in 
#£1122-1648 

The system at z = 2.030096 system in HE1 122-1648 shows 
C IV, N V, and possibly O VI absorption, while Lya is very 
weak (logAf(HI) = 1 1 .7 ± 0.3). This contrasts with most of the 
reliable systems we have found, where Lya is strong and N V is 
not detected. The line widths are such that a high temperature 
is possible, and so collisional ionization could be important. 

The density inferred from the photoionization models 
(log«H = -4.11 for model HM96, see Table f7|) is higher than 
one would expect for a photoionized, gravitationally confined 
cloud with the low H I column density seen. The neutral hydro- 
gen fraction log(nHi/«H) = -4.59, so with this density the cloud 
is very small: ~ 10 pc, and has a characteristic mass of only 
~1 Mq. Using equation (8) of Schaye (2001), we see that for 
the cloud to be self-gravitating, the gas fraction would have to 
be extremely small, f g ~ 10- 7 (n H /10" 411 cnr 3 r 3 r 2 2 :r 4 a52 , and 
the total mass would thus be of order 10 7 M Q . If the gas frac- 
tion is significantly greater than this, then an external pressure 
of P/k ~ 1 cm" 3 K would be required to confine the gas. This 
system may thus be a high redshift analogue of the low redshift 
system of Tripp et al. (2002) or the weak Mg II systems studied 
by Rigby, Charlton, & Churchill (2002). 

Since the ionization balance of a collisionally ionized plasma 
depends only on the temperature, we cannot constrain the den- 
sity if the system is collisionally ionized. Nevertheless, we 
can check whether the cloud can be self-gravitating by using 
the constraint that the differential Hubble flow across the ab- 
sorber must be smaller than the observed line width. Com- 
bining equation (ph with equation (4) of Schaye (2001), which 
gives nniNn, T, ff) for a self-gravitating cloud, yields the Hub- 
ble flow velocity across the system, 

For collisional ionization at a temperature T = 2.35 x 10 5 K 
the neutral hydrogen fraction is 10~ 5 7 , so for this system Nu = 
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10 17 - 4 cm" 2 . Thenb n - 1.2 x 10 4 /z kms _1 (/g/0.16). Thus, such 
a collisionally ionized system can only be self-gravitating if the 
gas fraction is extremely small. Note that if bn is much greater 
than the sound speed (~ 80 kms" 1 in this case), the absorber 
size exceeds the sound horizon, thermal pressure is irrelevant 
and the cloud cannot be in local, hydrostatic equilibrium (see 
Schaye2001). 

Both collisional and photoionization models (see Table ^ re- 
quire that the heavy element abundances be high - somewhere 
in the range 0.5 to 5 x solar. This system is also interesting in 
that the nitrogen abundance is high, and may even be enhanced 
relative to that of carbon and oxygen. 

In summary, regardless of whether the system is photoion- 
ized or collisionally ionized, its metallicity is high and it can 
only be self-gravitating if it has a negligible gas fraction. Alter- 
natively, it could be confined by external pressure or it could be 
a freely expanding transient. 

6. COSMOLOGICAL IMPLICATIONS 

An estimate the redshift evolution of the number of O VI 
systems between z ~ and z ~ 2 may be made by using the 
Tripp et al. (2000) low redshift results in conjunction with the 
data presented here. Tripp et al. find that the number per unit 
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48 for those systems in which the rest equiva- 



lent width W r > 30 mA for both OVI lines in the redshift range 
0.17 < z < 0.27. With 90% confidence they set > 17. The 
systems with redshifts z = 2.101764 and 2.352695 in HE1 122- 
1648, and 2.074740 and 2.075446 in HE22 17-28 18 satisfy their 
equivalent width criterion. This gives = 5.7 ±2.9, which is 
considerably lower than Tripp et al's value. The spectral reso- 
lution is comparable, but the S/N in the Tripp et al. data is not 
as high, so some of the blends separated in tables ^ and || might 
have been detected as single systems. Likely candidates are the 
close pairs at 2.080 and at z = 2.2145 in HE1 122-1648, and the 
complex at z=2. 180-2. 182 in 2217 would provide two systems 
above the threshold, as blends of the lowest two and highest 
two redshifted OVI's in the complex. Then there are a possible 
8 O VI systems above the threshold, so ^ = 1 1 .4 ± 4.0. This 
is consistent with the number density of such systems found in 
the redshift range 1.21 < z < 1.67 by Reimers et al. (2001), 
and strengthens the suggestion that the number density of O VI 
systems decreases with increasing redshift. However, there are 
large uncertainties in the number densities in the samples so 
further investigation is needed. 

It is interesting to check whether the identified O VI systems 
contain a cosmologically significant fraction of metals and/or 
baryons. The O VI density relative to the critical density is given 
by 
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For two redshift ranges centered on z = 2.2, Az = 0.365 and 
0.360, we have a mean total O VI column density per unit 
redshift EAfovi/Az = 1.26 x 10 15 , and so vi = 6 x 10 _8 //i- 
For the photoionization models log ^ ~ 0.75, so the cosmic 
metallicity at z = 2.2 due to O VI systems is 
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where we used Vt m = 0.3, D,\ = 0.7, and ili,h 2 = 0.02. The 
baryon density associated with the detected O VI systems is 
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The CLOUDY models give the total hydrogen/H I ratio 
«h/«hi ~ 10 5 . We do not have H I measurements directly as- 
sociated with the O VI component for two of the systems in 
HE2217-2818, so we base the estimate on the systems towards 
HE1 122- 1648. Then over a redshift range Az = 0.365 the total 
H I column density in systems for which O VI has been identi- 
fied is N(Hl) = 6.6 x 10 15 cm" 2 , and so at z = 2.2, 



^g.ovi ~ 6 x 10 3 h 



(7) 



For HE2217-2818 the total H I column density which may be 
associated with O VI systems is 1.9 x 10 16 cm" 2 . If this quantity 
is more appropriate, then £! g ,ovi is three times higher than the 
estimate above. Also, if the lower H I column density systems 
also contain O VI below the detection limit, then the value will 
be higher. In any case this is a significant fraction of the total 
baryo n density = 0.02/T 2 inferred from big bang nucleosyn- 



thesis O'Meara et al. 1997. Hence, we conclude that at z ~ 2 
the O VI systems contain a cosmologically significant fraction 
of both the metals and the baryons in the universe. 

7. CONCLUSIONS 

Specific O VI absorption systems in the 2 <; z <; 2.35 Lyman 
absorption forests of QSOs HE 1 122- 1648 and HE2217-2818 
were identified by self-consistently Voigt profile fitting the 
whole of the available spectral range (with all higher order Ly- 
man lines included). If these two QSO sight lines are typical, 
then at redshifts z < 2.5 O VI is usually present in absorption 
systems for which logV(HI) ;> 14, This is similar to the column 
density range for which C IV is detectable in spectra of z > 3 
QSOs (e.g., Songaila & Cowie 1996; Ellison et al. 2000). Weak 
O VI lines associated with lower column density Lya systems 
may also be present (Schaye et al. 2000), but are too difficult to 
detect directly. These observations indicate that O VI systems 
contain a significant fraction of the baryons (and metals) in the 
universe. 

For most systems the observed line widths, particularly those 
of the Lyman lines, do not allow temperatures high enough that 
O VI may be produced through collisional ionization. It is pos- 
sible to contrive a situation where the H I in collisionally ion- 
ized O VI regions is almost invariably blended with a cooler 
system with higher H I column density and no detected heavy 
elements, but even then collisional ionization would require that 
all systems have nearly the same unnatural temperature. 

For most systems photoionization models with a substantial 
hard ultraviolet flux reproduce the observations for expected 
densities and heavy element abundances generally in the range 
~ 10~ 3 -10~ 2 solar. These models give density estimates which 
are in good agreement with those from the photoionized, lo- 
cal hydrostatic equilibrium model of Schaye (2001). Assuming 
a hydrogen ionization rate Thi ~ 10~ 12 s" 1 at 2 <J z < 2.3, the 

hydrogen densities are generally within 0.5 dex of «h = 10 -4 . 

Photoionization by a UV spectrum much softer than that of 
Haardt & Madau (1996) (r H i/r H eii ~ 10 2 ) would yield lower 
densities than would be expected for self-gravitating clouds. 
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However, this possibility is ruled out because it would result 
in such large absorbers that the differential Hubble flow across 
the absorbers would far exceed the observed line widths. 

For some systems the interpretation is less straightforward. 
For example, a high ionization system towards HE1 122-1648 
at redshift z = 2.030 appears to have heavy element abundances 
close to solar, so is likely to be close to a metal-enriched galaxy. 
Intriguingly, regardless of whether the system is photoionized 
or collisionally ionized, it can only be self-gravitating if its gas 
fraction is negligible. Alternatively, it could be confined by 
external pressure or it could be freely expanding. In addition 
there are some absorption systems for which the C IV and O VI 
redshifts clearly do not agree well, and so there is some veloc- 
ity separation of the components in which these ions dominate. 
The system with the highest H I column density, at z = 2.180 



towards HE2217-2818, is the best example. Similar systems 
with high H I column density have been found by Simcoe et 
al. (2002). The relative frequency of occurrence of the various 
systems remains a topic for further investigation. 
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Table 1 

voigt profile parameters for systems in he 1122-1648 



Ion 


Redshift 


± 


b 


± 


logiV 


± 


HI 


2.100631 


000014 


20.84 


1.64 


13.392 


0.097 


HI 


2.101378 


000? 0? 


41.48 


11.19 


14.462 


0.351 


HI 


2.101764 


000006 


32.38 


4.36 


14.878 


0.172 


CIV 


2. 101764 




29.32 


1.69 


12.691 


0.017 


OVI 


2.101764 




29.25 


1.03 


14.127 


0.024 


HI 


2.102143 


00006 8 


24.64 


3.48 


14.030 


0.323 


HI 


2.102733 


0.000010 


25.58 


0.88 


13.452 


0.020 


HI 


2.104169 


0.000234 


32.19 


11.06 


12.748 


0.491 


HI 


2.104534 


0.000026 


29.18 


1.38 


13.554 


0.077 


HI 


2.106061 


0.000001 


19.98 


0.18 


13.676 


0.004 


HI 


2.106497 


0.000048 


151.55 


6.70 


13.368 


0.017 


HI 


2.107762 


0.000001 


21.74 


0.16 


13.922 


0.004 


HI 


2.109189 


0.000005 


15.96 


0.75 


12.685 


0.017 


HI 


2.110505 


0.000019 


26.02 


1.35 


13.239 


0.047 


HI 


2.110962 


0.000017 


29.62 


1.50 


13.409 


0.032 


HI 


2.112169 


0.000008 


35.69 


0.46 


13.921 


0.014 


HI 


2.112405 


0.000004 


15.57 


0.76 


13.505 


0.033 


HI 


2.113518 


0.000393 


39.75 


23.21 


12.879 


0.580 


HI 


2.114065 


0.000004 


29.59 


0.61 


14.932 


0.023 


HI 


2.115206 


0.000002 


32.58 


0.21 


14.879 


0.010 


HI 


2.116542 


0.000002 


19.77 


0.27 


13.320 


0.005 


HI 


2.117804 


0.000006 


22.81 


0.89 


12.797 


0.014 


HI 


2.120450 


0.000017 


25.41 


2.35 


12.398 


0.033 


HI 


2.121509 


0.000008 


36.72 


1.13 


12.972 


0.011 


HI 


2.125407 


0.000033 


75.49 


4.69 


12.689 


0.023 


HI 


2.128974 


0.000002 


20.01 


0.22 


13.433 


0.004 


HI 


2.130186 


0.000002 


18.39 


0.21 


13.436 


0.005 


HI 


2.130852 


0.000006 


32.07 


1.00 


13.143 


0.010 


HI 


2.131843 


0.000046 


33.06 


9.42 


12.229 


0.100 


HI 


2.132466 


0.000004 


21.43 


0.58 


13.245 


0.017 


HI 


2.133359 


0.000044 


52.87 


8.24 


13.147 


0.080 


HI 


2.133449 


0.000008 


24.10 


2.76 


12.992 


0.137 


HI 


2.136924 


0.000016 


27.50 


2.23 


12.587 


0.029 


HI 


2.142052 


0.000005 


25.05 


0.70 


12.928 


0.10 


HI 


2.144270 


0.000212 


55.27 


10.65 


13.174 


0.222 


HI 


2.144525 


0.000012 


23.53 


0.89 


14.065 


0.94 


HI 


2.144809 


0.000372 


30.76 


23.38 


13.124 


0.766 


HI 


2.145409 


0.000012 


20.41 


1.13 


12.995 


0.46 


HI 


2.146187 


0.000004 


27.39 


0.71 


13.129 


0.08 


HI 


2.146702 


0.000014 


13.99 


1.95 


12.168 


0.56 


HI 


2.149243 


0.000037 


49.08 


5.01 


12.428 


0.38 



Notes: Errors are la estimates. The Doppler parameter b is in km s" 1 and the ion column density N is the number per cm" 2 . Where 
no errors are given for the redshift its value is the same as for the ion above with the same quoted redshift value. The full table is 
available with the electronic version of this paper. 
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O VI at z ~ 2 



Table 2 

O VI SYSTEMS IN HE22 17-28 18 WITH 2.00 < z < 2.36 



Ion 


Redshift 


± 


b 


± 


o 


± 


^tuib 


T(K) 


HI 


2.074740 


0.000004 


28.07 


2.08 


14.890 


0.036 


14.2±0.9 


3.6±0.8x 10 4 


CIV 








15.87 


1.11 


12.887 


0.015 






OVI 








15.48 


0.69 


14.046 


0.155 






HI 


2.075398 


0.000003 


32.62 


6.90 


13.757 


0.472 


5.9±6.6 


6.3±3.0x 10 4 


CIV 








10.99 


0.51 


12.812 


0.015 






HI 


2.075446 


0.000006 


27.94 


3.83 


13.891 


0.290 


12.1 ±1.20 


3.9±1.4x 10 4 


OVI 


2.075446 




13.62 


0.87 


14.301 


0.029 






OVI 


2. 


180424 


0.000017 


3.71 


4.19 


12.781 


0.250 






CIV 


2. 


180488 


0.000149 


14.29 


10.21 


12.148 


0.603 






HI 


2. 


180546 


0.000382 


27.33 


9.06 


14.576 


1.296 






OVI 


2. 


180650 


0.000027 


25.81 


2.93 


13.723 


0.049 






CIV 


2. 


180739 


0.000026 


13.93 


1.61 


13.348 


0.093 






SilV 


2. 


180746 


0.000013 


12.78 


0.71 


12.697 


0.056 






Silll 


2. 


180747 


0.000011 


12.42 


0.65 


12.889 


0.054 






HI 


2. 


180791 


0.000011 


28.56 


7.32 


15.766 


0.397 






CII 


2. 


180791 




10.36 


1.68 


12.732 


0.050 






SilV 


2. 


180838 


0.000004 


4.82 


1.25 


12.315 


0.124 






Silll 


2. 


180843 


0.000006 


4.77 


1.34 


12.485 


0.108 






CIV 


2. 


180845 


0.000009 


7.00 


2.75 


12.682 


0.385 






CIV 


2. 


181155 


0.000031 


12.78 


4.44 


12.184 


0.122 






SilV 


2. 


181174 


0.000018 


3.85 


3.62 


11.336 


0.156 






HI 


2. 


181348 


0.000356 


49.59 


11.52 


14.841 


0.320 






OVI 


2. 


181371 


0.000022 


14.34 


3.29 


13.228 


0.079 






CIV 


2. 


181726 


0.000016 


8.87 


1.52 


12.450 


0.102 






OVI 


2. 


181781 


0.000015 


9.30 


2.14 


13.416 


0.134 






SilV 


2. 


181804 


0.000021 


6.60 


3.52 


11.455 


0.137 






CIV 


2. 


181917 


0.000075 


10.96 


8.50 


11.931 


0.347 






HI 


2. 


181978 


0.000060 


23.15 


4.51 


13.973 


0.263 






OVI 


2. 


182059 


0.000030 


19.67 


3.61 


13.723 


0.072 







Notes: The turbulent component of the Doppler parameter b tm \, and the temperature T are evaluated using the b 2 = b 2 mbujeBt + b 2 h 
model described in the text. Other details are as for Table 0. 
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Table 3 

O VI systems in HE1 122-1648 with 2.0 < z < 2.365 



Ion 


I 


± 


b 


± 


\ogN 


± 


blurb 


T(K) 


HI 


o nm i no 


fl AAAAAzl 


?9 46 


1 40 


1 3 .4 J 3 


099 


14 i 1 


A 1 -4- A « v 1 A 4 
4.1 I U.O X 1U 


CIV 






1 5 87 


47 


13 1 00 


014 






SilV 






14^82 


L13 


1 1.987 


o'024 






HI 


2.007205 


0.000005 


23.50 


2.49 


15.135 


0.302 


2±5 


3.4 ±0.8 x 10 4 


CIV 






6.98 


1.61 


12.309 


0.072 






OVI 






6.10 


1.60 


13.424 


0.080 






HI 


o a^iaaq^ 


n aaaaa^ 

U.UUUUUJ 


22 96 


5 1 2 


1 1 738 


300 


7 -1- 1 A 
/ I 1U 


a n -Hi s v i a 4 

U.U ZTL o.J X 1U 


CIV 






7 44 


96 


1 1 1 (r)^ 


040 






NV 






7 '44 


57 


12848 


021 






OVI 






7.44 


1.89 


13.385 


0.075 






HI 


2 033133 


n nnnnm 

U . UWUWU 1 


34 03 


25 


14 684 


011 


15 i 2 


< o _i_ a ri v i n 4 

J.O ^ X.J.J X 1 u 


OVI 






16 59 


1 4li 


1 3 678 


033 






HI 


2.064412 


0.000041 


39.00 


5.27 


14.241 


0.168 


39 ±6 


0.0 ±4.0 x 10 4 


CIV 






39.00 


5.56 


12.559 


0.070 






HI 


2.064789 


0.000016 


72.50 


21.46 


14.052 


0.420 


7 ± 17 


3.2 ±2.1 x 10 


CIV 






21.91 


2.19 


12.520 


0.070 






HI 


2.065681 


0.000014 


40.69 


11.78 


13.345 


0.266 


13 ± 13 


9.1 ±6.3 x 10 4 


CIV 






17.14 


5.58 


12.198 


0.040 






OVI 






16.21 


4.35 


13.029 


0.103 






HI 


i aqai at 


fl AAAA 1 Q 
U.UUUUlo 


zz. /u 


1 fY3 
1 AJJ 


t a 7nn 
14. /uu 


n nan 


OA -U 1 


A 7 ± 7 (1 v 1 A^ 
u./± i.UX 1U 


L 1 V 






on AA 


6.07 


IZAUZ 


0.057 






OVI 






20.38 


2.39 


13.503 


0.054 






HI 


2.080527 


0.000021 


17.74 


1.40 


14.236 


0.154 


IS ± 3 


0.0 ± 6.9 x 10 3 


CIV 






17.74 


8.39 


11.864 


0.093 






OVI 






17.74 


2.63 


13.354 


0.071 






HI 


2.101764 


0.000006 


32.38 


4.36 


14.878 


0.172 


29 ±1 


1.3 ± 1.9 x 10 4 


CIV 






29.32 


1.69 


12.691 


0.017 






OVI 






29.25 


1.03 


14.127 


0.024 






HI 


2.206460 


0.000005 


31.11 


0.54 


15.129 


0.012 


7±2 


5.6 ±0.3 x 10 4 


CIV 






11.23 


6.03 


11.948 


0.130 






OVI 






10.52 


1.35 


13.123 


0.097 






HI 


2.215271 


0.000005 


26.1 


0.S1 


14.711 


0.039 


3±5 


4.1 ±0.3 x 10 4 


CIV 






8.18 


1.89 


11.869 


0.059 






OVI 






7.27 


1.32 


13.077 


0.082 






HI 


2.215494 


0.000031 


38.40 


3.60 


14.335 


0.157 


38 ± 10 


0.3 ±4.9 x 10 4 


CIV 






37.78 


6.95 


12.113 


0.064 






OVI 






37.77 


9.34 


13.351 


0.164 






HI 


2.339223 


0.000154 


23.81 


2.90 


13.959 


0.943 






HI 


2.339260 


0.000007 


19.26 


10.51 


13.492 


2.446 


8±3 


1.9 ±2.6 x 10 4 


CIV 






9.48 


0.92 


11.932 


0.041 






OVI 






9.13 


1.51 


12.986 


0.047 






HI 


2.339386 


0.000283 


24.34 


9.21 


13.793 


1.060 






HI 


2.352170 


0.000005 


21.14 


0.78 


14.415 


0.038 


10 ± 1 


2.1 ±0.2 x 10 4 


CIV 


2.352170 




11.15 


0.79 


12.383 


0.024 






OVI 


2.352170 




10.81 


0.82 


12.770 


0.126 






CIV 


2.352464 


0.000007 


7.53 


1.13 


12.150 


0.058 






HI 


2.352695 


0.000015 


47.74 


3.17 


14.870 


0.249 


47 ±2 


0.5 ±2.3 x 10 4 


OVI 


2.352695 




46.97 


1.85 


13.971 


0.015 






HI 


2.352716 


0.000061 


57.30 


7.96 


14.549 


0.505 






CIV 


2.352762 


0.000005 


16.66 


0.86 


12.800 


0.020 






CIV 


2.353383 


0.000016 


33.96 


2.53 


12.667 


0.027 







Note: Details as for Table |[ 
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O VI at z ~ 2 



Table 4 

photoionization models for he1 122- 1648 systems 



Redsh ft 

z 


HI 


Log column density for: 
CIV NV OVI 


SilV 


IVInrlpl 

(see text) 


logfZH 


[C/H] 


Temperature(K) 


Z.UoUlOZ 


14.70 


12 10 




13 50 






-4.0 




6. / ± /.0 X 10 




14.70 


loin 
1Z. lu 


1Z.ZJ 


1 J.JU 


6.58 


HM96 


-4.15 


-2.28 


5.1 x 10 4 




14.70 


12.10 


12.27 


13.50 


7.55 


P1.5B0 


-3.83 


-2.14 


4.4 x 10 4 




14.70 


12.10 


12.28 


13.50 


7.20 


P1.5B 1 


-4.79 


-3.12 


6.0 x 10 4 




14.70 


12.10 


12.26 


13.50 


6.33 


P1.5B2 


-5.85 


-4.13 


4.7 x 10 4 




14.70 


12.10 


12.29 


13.50 


7.85 


P1.8B0 


-4.10 


-2.41 


4.7 x 10 4 




14.70 


12.10 


12.30 


13.50 


7.43 


P1.8B1 


-5.07 


-3.40 


5.9 x 10 4 




14.70 


12.10 


12.26 


13.50 


7.43 


P1.8B2 


-6.19 


-4.44 


4.1 x 10 4 




14.70 


12.07 


12.50 


13.53 


10.32 


Coll 




-2.89 


2.5 x 10 5 


2.080527 


14.24 


1 1 86 




1 3 35 






-4.3 




n a l /r n . , i r\i 

u.u ± o.y X 1U 




14.24 


1 1.86 


T7T — 

12.06 


13.35 


6.11 


HM% 


-4.22 


-2.00 


5.1 x 10 4 




14.24 


1 1 .50 


iz.u / 


n 

l j.jj 


7.07 


P1.5B0 


-3.91 


-1.89 


4.4 x 10 4 




14.24 


11.86 


12.08 


13.35 


6.77 


P1.5B 1 


-4.85 


-2.86 


6.1 x 10 4 




14.24 


11.86 


12.06 


13.35 


5.91 


P1.5B2 


-5.91 


-3.87 


4.6 x 10 4 




14 24 


11.86 


12.08 


13.35 


7 40 


PI 8B0 


-4 17 


-2 17 


4 8 x 10 4 




14.24 


11.86 


12.09 


13.35 


7.02 


P1.8B1 


-5.14 


-3.15 


6.0 x 10 4 




14.24 


11.86 


12.06 


13.35 


6.18 


P1.8B2 


-6.23 


-4.16 


4.0 x 10 4 




14.24 


11.87 


12.31 


13.34 


10.12 


Coll 




-2.63 


2.5 x 10 5 


2.101764 


14.88 


12.69 




14.13 






-3.9 




1.2 ± 1.9 x 10 4 




14.88 


12.69 


12.86 


14.13 


7.06 


HM96 


-4.19 


-1.85 


4.9 x 10" 




14.88 


12.69 


12.88 


14.13 


8.08 


P1.5B0 


-3.86 


-1.72 


4.2 x 10 4 




14.88 


12.69 


12.88 


14.13 


7.68 


P1.5B1 


-4.82 


-2.69 


5.9 x 10 4 




14.88 


12.69 


12.86 


14.13 


6.74 


P1.5B2 


-5.91 


-3.72 


4.8 x 10 4 




14.88 


12.69 


12.89 


14.13 


8.36 


P1.8B0 


-4.14 


-1.99 


4.6 x 10 4 




14.88 


12.69 


12.89 


14.13 


7.88 


P1.8B1 


-5.12 


-2.99 


5.8 x 10 4 




14.88 


12.69 


12.85 


14.13 


6.87 


P1.8B2 


-6.23 


-4.04 


4.2 x 10 4 




14.88 


12.69 


12.90 


14.13 


8.56 


P2B0 


-4.32 


-2.17 


4.8 x 10 4 




14.88 


12.69 


12.90 


14.13 


8.03 


P2B1 


-5.31 


-3.17 


5.8 x 10 4 




14.88 


12.68 


13.11 


14.14 


10.92 


Coll 




-2.47 


2.5 x 10 5 


2.206460 


15.13 


11.95 




13.12 






-3.7 




5.6 ±0.3 x 10 4 




15.13 


11.94 


11.99 


13.12 


6.89 


HM96 


-4.02 


-2.98 


4.9 x 10" 




15.13 


11.95 


12.03 


13.12 


7.94 


P1.5B0 


-3.49 


-2.78 


4.1 x 10 4 




15.13 


11.95 


12.04 


13.12 


7.42 


P1.5B1 


-4.65 


-3.77 


5.8 x 10 4 




15.13 


11.95 


12.00 


13.12 


6.43 


P1.5B2 


-5.73 


-4.82 


5.0 x 10 4 




15.13 


11.95 


12.06 


13.12 


8.16 


P1.8B0 


-3.94 


-3.03 


4.5 x 10 4 




15.13 


11.95 


12.06 


13.12 


7.61 


P1.8B1 


-4.94 


-4.03 


5.8 x 10 4 




15.13 


11.96 


12.40 


13.11 


11.18 


Coll 




-3.50 


2.35 x 10 5 




15.13 


11.95 


11.61 


13.12 


9.48 


OP1.8B0 


-3.39 


-3.03 


3.4 x 10 4 




15.13 


11.96 


11.68 


13.12 


8.79 


OP1.8B1 


-4.44 


-4.07 


5.0 x 10 4 



Notes: The quantities in the rows in which a redshift is given are the measured values, with the hydrogen density estimate «h (cm -3 ) 
and temperature derived as described in the text. The values below those lines give the results of CLOUDY models, as described in 
the text, in which the hydrogen density and heavy element abundances relative to solar ([C/H]) were varied to provide a best fit to the 
H I, C IV and O VI column densities. The temperatures there are the CLOUDY estimates. For collisional ionization models ("Coll" 
above) the temperature and abundances were varied to get the best fit to the three column densities. 
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Table 5 

voigt profile fits for collisional ionization models with imposed high temperatures 



Ion 


Redshift 


± 


b 


± 


logJV 


± 


Temperature (K) 


HI 


2 080268 


000002 


29.82 


0.26 


14.990 


0.016 




HI 


2 0801 52 


000032 


64.25 






13.70 


2.5 x 10 5 


CIV 


2 0801 52 




18.60 




11.684 


0.309 




OVI 


2.080152 




16.12 


7.10 


12.959 


0.355 




HI 


2.080384 


0.000047 


72.66 






13.26 


2.5 x 10 5 


CIV 


2.080384 




38.76 




12.240 


0.101 




OVI 


2.080384 




37.63 


2.93 


13.710 


0.070 




HI 


2.101668 


0.000005 


41.53 


1.09 


15.078 


0.014 




HI 


2.101769 


0.000008 


66.56 


0.51 


14.128 


0.132 


2.5 x 10 5 


CIV 


2.101769 




25.54 




12.666 


0.020 




OVI 


2.101769 




23.80 




13.875 


0.069 




HI 


2.206455 


0.000007 


29.26 


0.53 


15.144 


0.018 




HI 


2.206405 


0.000012 


62.27 




13.770 


0.169 


2.35 x 10 5 


CIV 


2.206405 




18.04 




12.451 


0.062 




OVI 


2.206405 




15.63 


1.48 


13.287 


0.178 





Notes: The temperature required to reproduce the C IV/O VI ratio was imposed on the H I at that redshift. The associated H I column 
density, and that of any new component required to fit the data are given. Other details as for Table 111. 



Table 6 

photoionization models for the z = 2.007 components in he1 122- 1648 



z 


HI 


CIV 


NV 


OVI 


SilV 


Model 


density 


[C/H] 


Temperature(K) 


2.007109 


15.43 


13.10 






11.99 




-3.5 




4.1 ±0.6 x 10 4 




15.43 
15.43 


13.10 
13.10 


11.54 
10.92 


11.13 
9.73 


11.99 
11.99 


P1.8B0 
P1.8B1 


-2.27 
-2.64 


-1.16 
-1.61 


1.9 x 10 4 
2.1 x 10 4 


2.007205 


15.14 


12.31 




13.42 






-3.7 




3.4±0.8 x 10 4 




15.14 
15.14 


12.31 
12.31 


12.38 
12.40 


13.42 
13.42 


8.63 
8.07 


P1.8B0 
P1.8B1 


-3.80 
-4.80 


-2.70 
-3.70 


4.3 x 10 4 
5.7 x 10 4 



Note: Description as for Table f| 



Table 7 

Models for the z = 2.030096 system in HE1 122- 1648 



- HI 


CIV 


NV 


OVI 


SilV 


Model 


l0g«H 


[C/H] 


Temperature(K) 


2.030096 11.74 


12.27 


12.85 


13.39 










0.0 ±8.5 x 10 4 


11.74 


12.27 


12.26 


13.39 


7.82 


HM96 


-4.11 


0.69 


2.0 x 10 4 


11.74 


12.27 


12.30 


13.39 


9.34 


P1.8B0 


-4.34 


0.70 


5.0 x 10 3 


11.74 


12.27 


12.33 


13.39 


8.64 


P1.8B1 


-4.96 


-0.30 


2.1 x 10 4 


11.74 


12.26 


12.71 


13.41 


10.49 


Coll 




0.20 


2.35 x 10 5 



Note: Description as for Table [| 



